MethCORR infers gene expression from DNA methylation and allows molecular analysis of ten common cancer types using fresh-frozen and formalin-fixed paraffin-embedded tumor samples.

BACKGROUND: Transcriptional analysis is widely used to study the molecular biology of cancer and hold great biomarker potential for clinical patient stratification. Yet, accurate transcriptional profiling requires RNA of a high quality, which often cannot be retrieved from formalin-fixed, paraffin-embedded (FFPE) tumor tissue that is routinely collected and archived in clinical departments. To overcome this roadblock to clinical testing, we previously developed MethCORR, a method that infers gene expression from DNA methylation data, which is robustly retrieved from FFPE tissue. MethCORR was originally developed for colorectal cancer and with this study, we aim to: (1) extend the MethCORR method to 10 additional cancer types and (2) to illustrate that the inferred gene expression is accurate and clinically informative. RESULTS: Regression models to infer gene expression information from DNA methylation were developed for ten common cancer types using matched RNA sequencing and DNA methylation profiles (HumanMethylation450 BeadChip) from The Cancer Genome Atlas Project. Robust and accurate gene expression profiles were inferred for all cancer types: on average, the expression of 11,000 genes was modeled with good accuracy and an intra-sample correlation of R2 = 0.90 between inferred and measured gene expression was observed. Molecular pathway analysis and transcriptional subtyping were performed for breast, prostate, and lung cancer samples to illustrate the general usability of the inferred gene expression profiles: overall, a high correlation of r = 0.96 (Pearson) in pathway enrichment scores and a 76% correspondence in molecular subtype calls were observed when using measured and inferred gene expression as input. Finally, inferred expression from FFPE tissue correlated better with RNA sequencing data from matched fresh-frozen tissue than did RNA sequencing data from FFPE tissue (P < 0.0001; Wilcoxon rank-sum test). CONCLUSIONS: In all cancers investigated, MethCORR enabled DNA methylation-based transcriptional analysis, thus enabling future analysis of cancer in situations where high-quality DNA, but not RNA, is available. Here, we provide the framework and resources for MethCORR modeling of ten common cancer types, thereby widely expanding the possibilities for transcriptional studies of archival FFPE material.

Publisher Correction: MethCORR modelling of methylomes from formalin-fixed paraffin-embedded tissue enables characterization and prognostication of colorectal cancer.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Enhanced Performance of DNA Methylation Markers by Simultaneous Measurement of Sense and Antisense DNA Strands after Cytosine Conversion.

BACKGROUND: Most existing DNA methylation-based methods for detection of circulating tumor DNA (ctDNA) are based on conversion of unmethylated cytosines to uracil. After conversion, the 2 DNA strands are no longer complementary; therefore, targeting only 1 DNA strand merely utilizes half of the available input DNA. We investigated whether the sensitivity of methylation-based ctDNA detection strategies could be increased by targeting both DNA strands after bisulfite conversion. METHODS: Dual-strand digital PCR assays were designed for the 3 colorectal cancer (CRC)-specific methylation markers KCNQ5, C9orf50, and CLIP4 and compared with previously reported single-strand assays. Performance was tested in tumor and leukocyte DNA, and the ability to detect ctDNA was investigated in plasma from 43 patients with CRC stages I to IV and 42 colonoscopy-confirmed healthy controls. RESULTS: Dual-strand assays quantified close to 100% of methylated control DNA input, whereas single-strand assays quantified approximately 50%. Furthermore, dual-strand assays showed a 2-fold increase in the number of methylated DNA copies detected when applied to DNA purified from tumor tissue and plasma from CRC patients. When the results of the 3 DNA methylation markers were combined into a ctDNA detection test and applied to plasma, the dual-strand assay format detected 86% of the cancers compared with 74% for the single-strand assay format. The specificity was 100% for both the dual- and single-strand test formats. CONCLUSION: Dual-strand assays enabled more sensitive detection of methylated ctDNA than single-strand assays.

MethCORR modelling of methylomes from formalin-fixed paraffin-embedded tissue enables characterization and prognostication of colorectal cancer.

Transcriptional characterization and classification has potential to resolve the inter-tumor heterogeneity of colorectal cancer and improve patient management. Yet, robust transcriptional profiling is difficult using formalin-fixed, paraffin-embedded (FFPE) samples, which complicates testing in clinical and archival material. We present MethCORR, an approach that allows uniform molecular characterization and classification of fresh-frozen and FFPE samples. MethCORR identifies genome-wide correlations between RNA expression and DNA methylation in fresh-frozen samples. This information is used to infer gene expression information in FFPE samples from their methylation profiles. MethCORR is here applied to methylation profiles from 877 fresh-frozen/FFPE samples and comparative analysis identifies the same two subtypes in four independent cohorts. Furthermore, subtype-specific prognostic biomarkers that better predicts relapse-free survival (HR = 2.66, 95%CI [1.67-4.22], P value < 0.001 (log-rank test)) than UICC tumor, node, metastasis (TNM) staging and microsatellite instability status are identified and validated using DNA methylation-specific PCR. The MethCORR approach is general, and may be similarly successful for other cancer types.

Characterization of genetic intratumor heterogeneity in colorectal cancer and matching patient-derived spheroid cultures.

Patient-derived in vitro cultures of colorectal cancer (CRC) may help guide treatment strategies prior to patient treatment. However, most previous studies have been performed on a single biopsy per tumor. The purpose of this study was to analyze multiple spatially distinct biopsies from CRCs and see how well intratumor heterogeneity (ITH) was recapitulated in matching patient-derived spheroids. Three to five biopsies were collected from six CRC tumors. Each biopsy was split in two; one half was used for spheroid culturing, while the other half was used for DNA and RNA purification. For two patients, lymph node metastases were analyzed. Somatic mutations were called from whole exome sequencing data. Each tumor contained mutations shared across all biopsies and spheroids, including major CRC drivers such as APC, KRAS, and TP53. At the same time, all tumors exhibited ITH on both mutation and copy number level. The concordance between biopsies and spheroids ranged between 40 and 70% for coding mutations. For three patients, the biopsy and spheroid from matching areas clustered together, meaning that the spheroid resembled the area of origin more than the other areas. However, all biopsies and spheroids contained private mutations. Therefore, multiple cultures from spatially distinct sites of the tumor increase the insight into the genetic profile of the entire tumor. Molecular subtypes were called from RNA sequencing data. When based on transcripts from both cancer and noncancerous cells, the subtypes were largely independent of sampling site. In contrast, subtyping based on cancer cell transcripts alone was dependent on sample site and genetic ITH. In conclusion, all examined CRC tumors showed genetic ITH. Spheroid cultures partly reflected this ITH, and having multiple cultures from distinct tumor sites improved the representation of the genetic tumor subclones. This should be taken into account when establishing patient-derived models for drug screening.

An Optimized Shotgun Strategy for the Rapid Generation of Comprehensive Human Proteomes.

This study investigates the challenge of comprehensively cataloging the complete human proteome from a single-cell type using mass spectrometry (MS)-based shotgun proteomics. We modify a classical two-dimensional high-resolution reversed-phase peptide fractionation scheme and optimize a protocol that provides sufficient peak capacity to saturate the sequencing speed of modern MS instruments. This strategy enables the deepest proteome of a human single-cell type to date, with the HeLa proteome sequenced to a depth of ∼584,000 unique peptide sequences and ∼14,200 protein isoforms (∼12,200 protein-coding genes). This depth is comparable with next-generation RNA sequencing and enables the identification of post-translational modifications, including ∼7,000 N-acetylation sites and ∼10,000 phosphorylation sites, without the need for enrichment. We further demonstrate the general applicability and clinical potential of this proteomics strategy by comprehensively quantifying global proteome expression in several different human cancer cell lines and patient tissue samples.

The non-coding variant rs1800734 enhances DCLK3 expression through long-range interaction and promotes colorectal cancer progression.

Genome-wide association studies have identified a great number of non-coding risk variants for colorectal cancer (CRC). To date, the majority of these variants have not been functionally studied. Identification of allele-specific transcription factor (TF) binding is of great importance to understand regulatory consequences of such variants. A recently developed proteome-wide analysis of disease-associated SNPs (PWAS) enables identification of TF-DNA interactions in an unbiased manner. Here we perform a large-scale PWAS study to comprehensively characterize TF-binding landscape that is associated with CRC, which identifies 731 allele-specific TF binding at 116 CRC risk loci. This screen identifies the A-allele of rs1800734 within the promoter region of MLH1 as perturbing the binding of TFAP4 and consequently increasing DCLK3 expression through a long-range interaction, which promotes cancer malignancy through enhancing expression of the genes related to epithelial-to-mesenchymal transition.

SNHG16 is regulated by the Wnt pathway in colorectal cancer and affects genes involved in lipid metabolism.

It is well established that lncRNAs are aberrantly expressed in cancer where they have been shown to act as oncogenes or tumor suppressors. RNA profiling of 314 colorectal adenomas/adenocarcinomas and 292 adjacent normal colon mucosa samples using RNA-sequencing demonstrated that the snoRNA host gene 16 (SNHG16) is significantly up-regulated in adenomas and all stages of CRC. SNHG16 expression was positively correlated to the expression of Wnt-regulated transcription factors, including ASCL2, ETS2, and c-Myc. In vitro abrogation of Wnt signaling in CRC cells reduced the expression of SNHG16 indicating that SNHG16 is regulated by the Wnt pathway. Silencing of SNHG16 resulted in reduced viability, increased apoptotic cell death and impaired cell migration. The SNHG16 silencing particularly affected expression of genes involved in lipid metabolism. A connection between SNHG16 and genes involved in lipid metabolism was also observed in clinical tumors. Argonaute CrossLinking and ImmunoPrecipitation (AGO-CLIP) demonstrated that SNHG16 heavily binds AGO and has 27 AGO/miRNA target sites along its length, indicating that SNHG16 may act as a competing endogenous RNA (ceRNA) "sponging" miRNAs off their cognate targets. Most interestingly, half of the miRNA families with high confidence targets on SNHG16 also target the 3'UTR of Stearoyl-CoA Desaturase (SCD). SCD is involved in lipid metabolism and is down-regulated upon SNHG16 silencing. In conclusion, up-regulation of SNHG16 is a frequent event in CRC, likely caused by deregulated Wnt signaling. In vitro analyses demonstrate that SNHG16 may play an oncogenic role in CRC and that it affects genes involved in lipid metabolism, possible through ceRNA related mechanisms.

Cellular disposal of miR23b by RAB27-dependent exosome release is linked to acquisition of metastatic properties.

Exosomes are small secreted vesicles that can transfer their content to recipient cells. In cancer, exosome secretion has been implicated in tumor growth and metastatic spread. In this study, we explored the possibility that exosomal pathways might discard tumor-suppressor miRNA that restricts metastatic progression. Secreted miRNA characterized from isogenic bladder carcinoma cell lines with differing metastatic potential were uncoupled from binding to target transcripts or the AGO2-miRISC complex. In metastatic cells, we observed a relative increase in secretion of miRNA with tumor-suppressor functions, including miR23b, miR224, and miR921. Ectopic expression of miR23b inhibited invasion, anoikis, angiogenesis, and pulmonary metastasis. Silencing of the exocytotic RAB family members RAB27A or RAB27B halted miR23b and miR921 secretion and reduced cellular invasion. Clinically, elevated levels of RAB27B expression were linked to poor prognosis in two independent cohorts of patients with bladder cancer. Moreover, highly exocytosed miRNA from metastatic cells, such as miR23b, were reduced in lymph node metastases compared with patient-matched primary tumors and were correlated with increments in miRNA-targeted RNA. Taken together, our results suggested that exosome-mediated secretion of tumor-suppressor miRNA is selected during tumor progression as a mechanism to coordinate activation of a metastatic cascade.

Putative cis-regulatory drivers in colorectal cancer.

The cis-regulatory effects responsible for cancer development have not been as extensively studied as the perturbations of the protein coding genome in tumorigenesis. To better characterize colorectal cancer (CRC) development we conducted an RNA-sequencing experiment of 103 matched tumour and normal colon mucosa samples from Danish CRC patients, 90 of which were germline-genotyped. By investigating allele-specific expression (ASE) we show that the germline genotypes remain important determinants of allelic gene expression in tumours. Using the changes in ASE in matched pairs of samples we discover 71 genes with excess of somatic cis-regulatory effects in CRC, suggesting a cancer driver role. We correlate genotypes and gene expression to identify expression quantitative trait loci (eQTLs) and find 1,693 and 948 eQTLs in normal samples and tumours, respectively. We estimate that 36% of the tumour eQTLs are exclusive to CRC and show that this specificity is partially driven by increased expression of specific transcription factors and changes in methylation patterns. We show that tumour-specific eQTLs are more enriched for low CRC genome-wide association study (GWAS) P values than shared eQTLs, which suggests that some of the GWAS variants are tumour specific regulatory variants. Importantly, tumour-specific eQTL genes also accumulate more somatic mutations when compared to the shared eQTL genes, raising the possibility that they constitute germline-derived cancer regulatory drivers. Collectively the integration of genome and the transcriptome reveals a substantial number of putative somatic and germline cis-regulatory cancer changes that may have a role in tumorigenesis.

Biological activity and biotechnological aspects of locked nucleic acids.

Locked nucleic acid (LNA) is one of the most promising new nucleic acid analogues that has been produced under the past two decades. In this chapter, we have tried to cover many of the different areas, where this molecule has been used to improve the function of synthetic oligonucleotides (ONs). The use of LNA in antisense ONs, including gapmers, splice-switching ONs, and siLNA, as well as antigene ONs, is reviewed. Pharmacokinetics as well as pharmacodynamics of LNA ONs and a description of selected compounds in, or close to, clinical testing are described. In addition, new LNA modifications and the adaptation of enzymes for LNA incorporation are reviewed. Such enzymes may become important for the development of stabilized LNA-containing aptamers.

Natural RNA circles function as efficient microRNA sponges.

MicroRNAs (miRNAs) are important post-transcriptional regulators of gene expression that act by direct base pairing to target sites within untranslated regions of messenger RNAs. Recently, miRNA activity has been shown to be affected by the presence of miRNA sponge transcripts, the so-called competing endogenous RNA in humans and target mimicry in plants. We previously identified a highly expressed circular RNA (circRNA) in human and mouse brain. Here we show that this circRNA acts as a miR-7 sponge; we term this circular transcript ciRS-7 (circular RNA sponge for miR-7). ciRS-7 contains more than 70 selectively conserved miRNA target sites, and it is highly and widely associated with Argonaute (AGO) proteins in a miR-7-dependent manner. Although the circRNA is completely resistant to miRNA-mediated target destabilization, it strongly suppresses miR-7 activity, resulting in increased levels of miR-7 targets. In the mouse brain, we observe overlapping co-expression of ciRS-7 and miR-7, particularly in neocortical and hippocampal neurons, suggesting a high degree of endogenous interaction. We further show that the testis-specific circRNA, sex-determining region Y (Sry), serves as a miR-138 sponge, suggesting that miRNA sponge effects achieved by circRNA formation are a general phenomenon. This study serves as the first, to our knowledge, functional analysis of a naturally expressed circRNA.

Engineering small interfering RNAs by strategic chemical modification.

Synthetic small interfering RNAs (siRNAs) have revolutionized functional genomics in mammalian cell cultures due to their reliability, efficiency, and ease of use. This success, however, has not fully translated into siRNA applications in vivo and in siRNA therapeutics where initial optimism has been dampened by a lack of efficient delivery strategies and reports of siRNA off-target effects and immunogenicity. Encouragingly, most aspects of siRNA behavior can be addressed by careful engineering of siRNAs incorporating beneficial chemical modifications into discrete nucleotide positions during siRNA synthesis. Here, we review the literature (Subheadings 1 -3) and provide a quick guide (Subheading 4) to how the performance of siRNA can be improved by chemical modification to suit specific applications in vitro and in vivo.

In vivo knockdown of Brachyury results in skeletal defects and urorectal malformations resembling caudal regression syndrome.

The T-box transcription factor BRACHYURY (T) is a key regulator of mesoderm formation during early development. Complete loss of T has been shown to lead to embryonic lethality around E10.0. Here we characterize an inducible miRNA-based in vivo knockdown mouse model of T, termed KD3-T, which exhibits a hypomorphic phenotype. KD3-T embryos display axial skeletal defects caused by apoptosis of paraxial mesoderm, which is accompanied by urorectal malformations resembling the murine uro-recto-caudal syndrome and human caudal regression syndrome phenotypes. We show that there is a reduction of T in the notochord of KD3-T embryos which results in impaired notochord differentiation and its subsequent loss, whereas levels of T in the tailbud are sufficient for axis extension and patterning. Furthermore, the notochord in KD3-T embryos adopts a neural character and loses its ability to act as a signaling center. Since KD3-T animals survive until birth, they are useful for examining later roles for T in the development of urorectal tissues.

Spatio-temporal regulation of ADAR editing during development in porcine neural tissues.

Editing by ADAR enzymes is essential for mammalian life. Still, knowledge of the spatio-temporal editing patterns in mammals is limited. By use of 454 amplicon sequencing we examined the editing status of 12 regionally extracted mRNAs from porcine developing brain encompassing a total of 64 putative ADAR editing sites. In total 24 brain tissues, dissected from up to five regions from embryonic gestation day 23, 42, 60, 80, 100 and 115, were examined for editing.   Generally, editing increased during embryonic development concomitantly with an increase in ADAR2 mRNA level. Notably, the Gria2 (GluR-B) Q/R site, reported to be ~100% edited in previous studies, is only 54% edited at embryonic day 23. Transcripts with multiple editing sites in close proximity to each other exhibit coupled editing and an extraordinary incidence of long-range coupling of editing events more than 32 kb apart is observed for the kainate glutamate receptor 2 transcript, Grik2. Our study reveals complex spatio-temporal ADAR editing patterns of coordinated editing events that may play important roles in the development of the mammalian brain.

Development of Therapeutic-Grade Small Interfering RNAs by Chemical Engineering.

Recent successes in clinical trials have provided important proof of concept that small interfering RNAs (siRNAs) indeed constitute a new promising class of therapeutics. Although great efforts are still needed to ensure efficient means of delivery in vivo, the siRNA molecule itself has been successfully engineered by chemical modification to meet initial challenges regarding specificity, stability, and immunogenicity. To date, a great wealth of siRNA architectures and types of chemical modification are available for promoting safe siRNA-mediated gene silencing in vivo and, consequently, the choice of design and modification types can be challenging to individual experimenters. Here we review the literature and devise how to improve siRNA performance by structural design and specific chemical modification to ensure potent and specific gene silencing without unwarranted side-effects and hereby complement the ongoing efforts to improve cell targeting and delivery by other carrier molecules.

MicroRNA alterations and associated aberrant DNA methylation patterns across multiple sample types in oral squamous cell carcinoma.

BACKGROUND: MicroRNA (miRNA) expression is broadly altered in cancer, but few studies have investigated miRNA deregulation in oral squamous cell carcinoma (OSCC). Epigenetic mechanisms are involved in the regulation of >30 miRNA genes in a range of tissues, and we aimed to investigate this further in OSCC. METHODS: TaqMan® qRT-PCR arrays and individual assays were used to profile miRNA expression in a panel of 25 tumors with matched adjacent tissues from patients with OSCC, and 8 control paired oral stroma and epithelium from healthy volunteers. Associated DNA methylation changes of candidate epigenetically deregulated miRNA genes were measured in the same samples using the MassArray® mass spectrometry platform. MiRNA expression and DNA methylation changes were also investigated in FACS sorted CD44(high) oral cancer stem cells from primary tumor samples (CSCs), and in oral rinse and saliva from 15 OSCC patients and 7 healthy volunteers. RESULTS: MiRNA expression patterns were consistent in healthy oral epithelium and stroma, but broadly altered in both tumor and adjacent tissue from OSCC patients. MiR-375 is repressed and miR-127 activated in OSCC, and we confirm previous reports of miR-137 hypermethylation in oral cancer. The miR-200 s/miR-205 were epigenetically activated in tumors vs normal tissues, but repressed in the absence of DNA hypermethylation specifically in CD44(high) oral CSCs. Aberrant miR-375 and miR-200a expression and miR-200c-141 methylation could be detected in and distinguish OSCC patient oral rinse and saliva from healthy volunteers, suggesting a potential clinical application for OSCC specific miRNA signatures in oral fluids. CONCLUSIONS: MiRNA expression and DNA methylation changes are a common event in OSCC, and we suggest miR-375, miR-127, miR-137, the miR-200 family and miR-205 as promising candidates for future investigations. Although overall activated in OSCC, miR-200/miR-205 suppression in oral CSCs indicate that cell specific silencing of these miRNAs may drive tumor expansion and progression.

miRNA-dependent gene silencing involving Ago2-mediated cleavage of a circular antisense RNA.

MicroRNAs (miRNAs) are ∼22 nt non-coding RNAs that typically bind to the 3' UTR of target mRNAs in the cytoplasm, resulting in mRNA destabilization and translational repression. Here, we report that miRNAs can also regulate gene expression by targeting non-coding antisense transcripts in human cells. Specifically, we show that miR-671 directs cleavage of a circular antisense transcript of the Cerebellar Degeneration-Related protein 1 (CDR1) locus in an Ago2-slicer-dependent manner. The resulting downregulation of circular antisense has a concomitant decrease in CDR1 mRNA levels, independently of heterochromatin formation. This study provides the first evidence for non-coding antisense transcripts as functional miRNA targets, and a novel regulatory mechanism involving a positive correlation between mRNA and antisense circular RNA levels.

Enhancing miRNA annotation confidence in miRBase by continuous cross dataset analysis.

The immaculate annotation of all microRNAs (miRNAs) is a prerequisite to study their biological function on a genome-wide scale. However, the original criteria for proper miRNA annotation seem unsuited for the automated analysis of the immense number of small RNA reads available in next generation sequencing (NGS) datasets. Here we analyze the confidence of past miRNA annotation in miRBase by cross-analyzing publicly available NGS datasets using strengthened annotation requirements. Our analysis highlights that a large number of annotated human miRNAs in miRBase seems to require more experimental validation to be confidently annotated. Notably, our dataset analysis also identified almost 300 currently non-annotated miRNA*s and 28 novel miRNAs. These observations hereby greatly increase the confidence of past miRNA annotation in miRBase but also illustrate the usefulness of continuous re-evaluating NGS datasets in the identification of novel miRNAs.

Chemical modification of small interfering RNA.

Chemically synthesized siRNAs are widely used for gene silencing. For in vitro applications, stability, delivery, and immunological issues are rarely problematic, but for in vivo applications the situation is different. Limited stability, undesirable pharmacokinetic behaviour, and unanticipated side effects from the immune system call for more careful structural siRNA design and inclusion of chemical modifications at selected positions. Also the notion that siRNA induces significant off-target silencing of many non-related genes has promted new effective measures to enhance specificity. The scope of this review is to provide a simple guide to successful chemical and structural modification of siRNAs with improved activity, stability, specificity, and low toxicity.